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Abstract 

A model for spin-charge separated superconductivity in two dimensions is introduced 
where the phases of the spinon and holon order parameters couple gauge-invariantly to 
a statistical gauge-field representing chiral spin-fluctuations. The model is analyzed in 
the continuum limit and in the low-temperature limit. In both cases we find that phys- 
ical electronic phase correlations show a superconducting-normal phase transition of the 
Berezinskii-Kosterlitz-Thouless type, while statistical gauge-field excitations are found to 
be strictly gapless. It is argued that the former transition is in the same universality class 
as that of the XY model. We thus predict a universal jump in the superfluid density at this 
transition. The normal-to-superconductor phase boundary for this model is also obtained 
as a function of carrier density, where we find that its shape compares favorably with that 
of the experimentally observed phase diagram for the oxide superconductors. 
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I. Introduction 



Although the phenomenon of high-temperature superconductivity 1 has sparked a 
tremendous amount of theoretical activity, a viable theory remains to be uncovered. It 
is, however, generally agreed upon that strong electron-electron interactions must play an 
important role in the charge dynamics of the Copper-Oxygen planes that are common to 
these systems. 2 Among the various attacks on the theoretical problem of strongly inter- 
acting electrons in two dimension, those based on the hypothesis of spin-charge separation 
show promise at the phenomenological level. 2-3 In particular, gauge theories for the un- 
conventional metallic states found in the t — J model, which is considered to be the simplest 
model containing the essential strong correlation physics of the oxide superconductors, 4 
successfully account for many of the unconventional transport and collective mode prop- 
erties that are common to the corresponding metallic phases of these materials. 5 ' 6 For 
example, the T-linear in-plane resistivity, as well as the paradoxical observations of a hole- 
type Hall effect in conjunction with a large Luttinger Fermi surface, can be accounted for 
by such theories. 

In a similar spirit, analogous Ginzburg-Landau theories of spin-charge separated su- 
perconductivity itself have recently been proposed in the literature. 7 ' 8 It is presumed in 
such theories that separate superfluid instabilities exists in both the spinon and holon sys- 
tems in isolation. Here, the spinon superfluidity is assumed to arise from singlet Cooper 
pairing d la the mean-field resonating valence bond (RVB) scenario, 9 while the holon su- 
perfluidity is driven by Bose-Einstein condensation. Given the Ioffe-Larkin composition 
law, R = Rb + Rf : which states that the physical resistance is given by the sum of the 
holon (6) and spinon (/) contributions, then true superconductivity occurs only when both 
species are superfluid. 3 However, because such Ginzburg-Landau theories are formulated 
in the continuum, the analysis of the superconducting transition in two dimensions is 
complicated by the structure of the vortex cores. In addition, these theories have been 
only analyzed in the mean-field approximation, to date, where fluctuations of the statisti- 
cal gauge-field related to the constraint against double occuppancy are neglected. (Note 
that such excitations physically represent so-called chiral spin fluctuations. 10 ) As a result, 
the precise shape that they obtain for the phase boundary between the normal and the 
superconducting state is questionable. 8 
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In order to address the role played by fluctuations of the statistical gauge-field, we 
introduce here an alternate spin-charge separated model of strongly correlated supercon- 
ductivity in which only the phase of each order parameter is allowed to vary. This model 
is a two-component generalization of a lattice gauge theory model in two euclidean (1+1) 
dimensions known as the Abelian-Higgs model. 11-13 Specifically, the theory contains a 
holon-pair phase order parameter (Higgs field) and a spinon-pair phase order parameter 
(Higgs field) that couple in a gauge-invariant way to the fluctuating U(l) (Abelian) sta- 
tistical gauge-field describing chiral spin fluctuations. 10 We henceforth refer to the latter 
as the two-component Abelian-Higgs (AH 2 ) model. 14 By considering both the contin- 
uum and the low temperature limits of this model, we find that it exhibits a normal-to- 
superconducting transition of the Berezinskii-Kosterlitz-Thouless (BKT) type in strictly 
two spacial dimensions. 15 ' 16 Note that the low-temperature analysis is achieved by a Vil- 
lain duality transformation of this model, 17 which has been successfully employed in the 
past to study the one-component Abelian Higgs model and the XY model. 11 ' 18 ' 19 We ob- 
tain the following two major results: that (i) the phase correlations corresponding to the 
statistical gauge- field are short-range at all temperatures, and that (ii) only phase correla- 
tions corresponding to the physical electronic order parameter show algebraic long-range 
order below a BKT- type transition temperature, T c = §(</ 6 _1 + J7" 1 ) -1 , where Jb and 
Jf are the respective local phase-stiffnesses of the holon-pairs and the spinon-pairs. The 
former result (i) implies that statistical gauge-field fluctuations do not acquire a gap via 
the Higgs mechanism. The latter result (ii), on the other hand, implies that this spin- 
charge separated model shows a true superconducting transition at T c . 16 It is important 
to remark that the present calculation, which incorporates fluctuations of the statistical 
gauge-field, results in a transition temperature inferior to the mean-field approximation 
result, 3 ' 8 Tc°^ = ^ min( J/), as expected. In addition, the presently obtained transition 
temperature yields a metal-superconductor phase diagram as a function of hole doping 
that qualitatively resembles that of the oxide superconductors (see Fig. I). 20 It is also 
argued that this transition falls into the same universality class as that of the XY model. 
We thus predict that the present AH 2 model for spin-charge separated superconductivity 
has a universal jump in the superfluid density at the transition. 16 Last, we compute the 
Wilson loop and find that it shows a perimeter law in the superconducting phase. On 
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the other hand, by continuity with the corresponding results obtained for the case of the 
one-component model, 11 we conclude that the Wilson loop generally shows a "confining" 
area law in the normal phase. This change in behavior simply reflects the vortex binding- 
unbinding transition. We then argue, however, that the latter "confinement" effect is 
a trivial result of "electromagnetism" in one space and one time dimension (1+1), and 
hence that fluctuations in the statistical gauge-field remain gapless in the normal phase. 21 
This, hence, provides a basis for the calculations of linear-in-T resistance in the "strange" 
metallic phase of the t — J model, 5 ' 6 which rely on this property. Note that gapless sta- 
tistical gauge-field excitations have also been shown to exist in a two-component anyon 
superconductor saddle-point of the t — J model known as the commensurate flux-phase 
at low temperature and near half-filling. 22 ' 23 Yet once dynamical effects are included, an 
exponentially small gap appears close to half-filling in such case. 22 

The remainder of the paper is organized as follows. In the next section we introduce 
the AH 2 model in the context of strongly interacting electrons in two dimensions. Section 
III contains a discussion of various limits of the model, including the continuum limit. 
In section IV we study the phase correlations in the low-temperatue limit via the Villain 
duality transformation, while the nature of the statistical gauge-field excitations are treated 
in section V within the same context. A renormalization group analysis of the BKT- 
transition found in this model is presented in section VI, while we compare the theoretically 
obtained phase diagram with that observed experimentally in the oxide superconductors 
as function of hole doping in section VII. The nature of the normal state of the present 
model is also discussed here. Last, the phase correlators of the model are computed in the 
"spin-wave" approximations in Appendix A, and those of the one-component model are 
computed in the low-temperature limit in Appendix B. 



II. Two-component Abelian-Higgs Model 

Before introducing the Abelian-Higgs model for spin-charge separated superconduc- 
tivity in two dimensions, let us first consider the corresponding normal state, which we 
take to be the unconventional metallic state found in strongly interacting two-dimensional 
(2D) electron systems first discussed by Ioffe and Larkin in terms of the spin-charge sepa- 
rated language of spinons and holons. 2 ' 3 The archetypical theoretical description of strong 
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correlation physics in the context of the oxide superconductors is given by the t — J model 
defined over the square lattice, 4 notably where double occuppancy of electrons at each site 
is excluded as a result such correlations. This constraint may be imposed by the intro- 
duction of an auxilliary slave-boson field, bi, such that the electron field is re-expressed 
as 

Q<7 * Q<7^j ? (1) 

along with the constraint c\ a Ci a + b\bi = 1 at each site. 22 Certain mean-field treatments of 
the latter constraint allowing the slave-boson to propagate result in a spin-charge separated 
description of the t — J model, with Cj CT representing the spinon field and bi representing 
the holon field. 3-6 A statistical gauge field, A^, generated by the internal gauge symmetry 
{cicr.bi) — > (e^Cjo-, e ldi bi) shown in Eq. (1) then mediates interactions between the two 
species. Within this representation, the "strange" metal is then characterized by a normal 
spinon Fermi-liquid state with a large Luttinger Fermi surface in conjunction with a normal 
(uncondensed) holon liquid state. The statistical gauge-field mediates interactions among 
both species. 5 ' 6 In particular, interactions of these fluctuation with the uncondensed holon 
liquid ultimately give rise to the T-linear prediction for the resistance characteristic of the 
"strange" metal phase. 5 

Upon lowering temperature, however, meanfield resonanting- valence-bond (RVB) stud-| 
ies of the t — J model show that amplitude develops for an order parameter describing a 
paired spinon state. 9 In addition, amplitude for a separate order parameter representing 
the Bose-Einstein holon condensate exists at low-temperature. Hence by definition, in the 
mean-field approximation, the unconventional metallic state can only exist at temperatures 
above both the highest temperature, Tf, at which any spinon Cooper pairing instability 
sets in, and above the Bose-Einstein condensation temperature, T b , for the ho Ions (see Fig. 
1). Therefore, as temperature, T, is lowered, the following three situations may arise: (i) 
Tf < T < Tb] (ii) Tb < T < Tf, and (Hi) T < Tb,Tf. Given the Ioffe-Larkin composition 
law, 3 R = Rb + Rf, which states that the physical resistance is given by the sum of the 
holon (b) and spinon (/) contributions, then only regime (in) is superconducting. This 
implies that the critical temperature for the superfluid transition is given by 3 

Ti°)=min(T b ,T / ) (2) 

at the mean-field level. As has been discussed in the literature, 8 regimes (i) and (ii) are 
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non-superconducting phases that correspond, respectively, to a normal Fermi-liquid state 
and a to spin-gap state. The various regimes mentioned above are shown in Fig. 1. 

The latter phase diagram for strongly interacting electrons in two dimensions is ob- 
tained in the mean-field approximation, however, where fluctuations of the statistical 
gauge-field are absent. In this paper, we shall study what effect such fluctuations have 
on the previously discussed superfluid transitions occurring in the spin-charge separated 
metallic state by focusing on the phase degrees of freedom alone of both the spinon and the 
holon superfluid order parameters. In particular, consider the following partition function 
for the two-component square lattice Abelian-Higgs model, 14 which we presume models 
the superfluid sector of the spinon/holon system in the presence of statistical gauge- field 
fluctuations: 

Z = J VMr)VM r ') V M r ">M-E/T), (3) 

where 

| =P b £{1 - cos[A M </> b (r) " qA^r)]} + (3 f £{1 - cos^^r) - qA^r)}}+ 

+ 2^2 E V ~ cos i^Mr) ~ A„4,(r)]}. ( 4 ) 

Here, 06 (r) and 4>f(r) represent the respective phases of the holon and spinon order param- 
eters, A^{r) = Ajtf+p, denotes the statistical gauge-field, 3-8 and A M denotes the lattice 
difference operator 11 (/z = x,y). The first two terms above are the lattice versions of the 
stiffness energies for the phase fluctuations, while the last term is the corresponding stiff- 
ness energy for the statistical gauge-field fluctuations. Also, (3 b = J b /T and (3f = Jf/T, 
where J b and J/ denote the local phase stiffness of the holon order parameter and the 
spinon order parameter, respectively, while g~ 2 = Xd/T, where Xd represents the local 
stiffness for statistical gauge-field fluctuations. (Throughout, we set ks, and the lattice 
constant, a, to unity, and we sum over repeated indices.) Spin-charge separated treatments 
of the t — J model in two dimensions give a diamagnetic susceptibility for fluctuations in 
the statistical gauge-field of Xd = Xf + Xb, where Xf an d Xb are the diamagnetic suscepti- 
bilities of each species. 5 ' 6 Within the metallic (fluxless) saddle-point of these treatments, 
Xf ~ J(l — x) and Xb ~ tT b /T, where x denotes the hole concentration and T b ~ tx denotes 
the ideal Bose-Einstein condensation temperature. Also, J b ~ fx, where t' ~ t denotes 
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the relevant matrix element for hopping in the holon liquid, and Jf ~ J(l — x/xo), where 
xo ~ 1 denotes the critical hole concentration above which the spinon pairing instability is 
absent. 9 Finally, q represents the charge per site of each species. (Throughout, we follow 
the notations used in refs. 11, 18, and 19.) We remind the reader that Landau- Ginzburg 
versions of this model, formulated in the 2D continuum, have already been discussed in 
the literature. 7 ' 8 

Before we begin the analysis of the AH 2 model, a few remarks are in order. First, since 
physical electrons have no statistical charge by Eq. (1), we have restricted our analysis 
(4) to statistically charge neutral AH 2 models; i.e., the statistical charge per site, q, of 
the holons and of the spinons is equal. Second, we implicitly presume that the spinon 
order parameter, <f>f, results from Cooper pairing of the spinons a la the RVB picture for 
superconductivity. 9 Hence, from this point on we shall take a charge per site of q = 2 for 
the present model, which implies that the holon condensate is made up pairs of holons 
as well. Note that a holon pairing instability consistent with the previous assignment has 
been shown to exist in the "strange" metallic phase of the t — J model, 6 which is the 
normal state of the present model. Also, in the limit near half-filling relevant to the oxide 
superconductors, the density of holon pairs is small and the overlap between two pairs 
of bosons should be negligible. Therefore, a bose condensate of dilute bosonic molecules 
can form, 24 resulting in the holon order parameter, above. If, on the other hand, the 
holon condensate where to result from a condensation of unpaired holons as is suggested 
by certain mean-field treatments, 4 ' 5 ' 9 then such an order parameter must necessarily be 
defined on a lattice different from that of the spinon order parameter by the requirement 
of statistical charge neutrality (1). The latter physical situation, however, cannot be 
described by the present AH 2 model, which has only one underlying lattice. 



III. Limits 

The nature of the above AH 2 model (4) for spin-charge separated superconductivity 
on the square lattice can be uncovered, in large part, by an analysis of various limits, as 
we discuss below. 

A. Continuum Limit. Consider the limit where q 2 g 2 (3b, q 2 g 2 (3f *C 1 and where fluctu- 
ations in the fields are at longwavelength, which corresponds to the continuum limit [see 
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Eqs. (4) and (20)]. In this case, the last term in the energy functional (4) is negligible, 
whereas the first two terms may approximated by 

% = \\ d 2 rMVcj> b - qAf + PfiVcPf - qA) 2 ]. (5) 

As is customary in the treatment of the XY model, let us now separate the phase con- 
figurations into "spin- wave" and vortex components; 19 i.e., let 4> a = (j)^ + (p^ for each 
species s = b, f, where 4>s represents the smooth portion of the configuration such that 
V x (V(pi w ^) = 0, and where <fi^ represents the portion of the configuration containing 
vortex singularities such that V • (V^i^) = 0. Taking the Coulomb gauge, V • A = 0, we 
see that statistical gauge-field fluctuations only communicate with the component of the 
phase configuration containing vorticity, since the statistical gauge-field is purely tranverse 
in this gauge, as are such configurations. Therefore, integrating out first the gauge-field 
excitations in the corresponding partition function (3) within this gauge leaves us with 
following effective action for the vortex component: 

^ = f/^ 2 r(V^-V^) 2 , (6) 

where 

^(P^+PJ 1 )- 1 - (7) 
The remaining smooth "spin-wave" components result in a trivial Gaussian action, 

Hence, by eqs. (6) and (7) we expect a BKT transition to occur in the physical electronic 
phase, 

0el = 0/ -06 (8) 
when 2/T/3 = 4, which implies a corresponding BKT transition temperature of 

Tc^l^ + Jj 1 )- 1 . (9) 

In the continuum limit, therefore, statistical gauge-field fluctuations suppress, but do not 
destroy, the superconductivity in this spin-charge separated system in relation to the Ioffe- 
Larkin mean- field approximation result; i.e. < T c < min(T;,,T/), where T b = and 
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Tf = ^Jf are the corresponding transition temperatures of each specie separately. As we 
will show in the next section, this result is generally valid outside of the continuum limit, 
as well. 

B. XY model + Pure Gauge Theory. Consider now the limit where (3b (or (3f) ^> 1. 
After making the gauge-transformation A'^ = — q~ 1 A^(j)b, we can rewrite the energy 
functional (4) as 

| =(3 b £[1 - cos(^y] + (3 f £{1 - cos[A M el (r) - qA'^r)]}+ 

r,n r,n 

+ E i 1 - C ^ A » - A„4.(r)]}. (10) 

Hence, the present limit imposes the constraint A'^ = 2irn/q, where n is any integer. The 
above energy functional then reduces to 

£{1 " cos[A M el (r)]} + ^{1- cos[A^(r) - A^(r)]}. (11) 

The latter expression, in conjunction with the former constraint, contains an XY model 
that is decoupled from a pure gauge theory known in the literature as Z q . 13 The XY 
model part implies a BKT transition at Ty for the physical electronic phase, e i, which is 
consistent with the results (9) obtained previously in the continuum limit for the present 
case, Jb — > oo. The (/-state (Z q ) pure gauge theory part of the energy functional (11), on 
the other hand, can be easily shown to exhibit no phase transition in two dimensions by 
choosing to work in the Landau gauge (A x = 0). 13 This, of course, reduces the theory to 
one dimension. 

C. One- component Abelian-Higgs Model. Last, let us take the limit (3b = or (3f = 0, 
which corresponds to the well studied one-component Abelian-Higgs model. The latter 
model sustains no Higgs phenomenon, since the long-range interaction is screened, and 
in general, it is thought to be free of phase transitions in two dimensions. 11-13 In fact, 
the phase correlation length can be shown to be finite in the low-temperature limit (see 
Appendix B). Again, these results are consistent with those obtained in the continuum 
limit (9); i.e., T c = in the case that Jb = or Jf = 0. 

In summary, we find that the physical electronic phase (8) undergoes a BKT tran- 
sition at a critical temperature given by expression (9), while the gauge-field excitations 
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experience no phase transitions at all, in the limits considered above. Below, we show that 
these findings are corroborated by a low-temperature analysis of the model. 

IV. Phase Correlations 

With the intent of understanding the superconducting properties of our model (4) 
for spin-charge separated superconductivity in two dimensions, let us consider now the 
low-temperature limit f3j,,f3f,g~ 2 ^> 1. In this limit, the most important configurations of 
both the phases and the statistical gauge-field lie at extrema of the cosine functions found 
in expression (4). Thus, we may employ well-known Villain duality transformations that 
are successful in analysing both the XY model 18 ' 19 and the one-component Abelian-Higgs 
model 11 in the same limit. This approximation amounts to replacing the exponential of a 
cosine function by the sum of exponentials of parabolas, each situated at the appropriate 
extrema of the cosine function. 17 Below, we show that the latter analysis leads to a gen- 
eralized Coulomb gas ensemble that exhibits a BKT phase transition consistent with the 
discussion found in the previous section (III. A). Note also that the analogous treatment 
of the phase auto-correlations function in the special case of the one-component Abelian- 
Higgs model ((3b = or j3f = 0), which is useful in the study of the unconventional normal 
state corresponding to the present model for spin-charge separated superconductivity, 21 is 
found in Appendix B. 

A. Coulomb Gas Representation. Application of the above mentioned low-temperature 
approximation to the presently considered AH 2 model reduces to the substitution of the 
mathematical identity 



that is valid in the limit j3 — > oo, into the partition function (3) corresponding to the energy 
functional (4). By closely following the analogous treat of the one-component version of 
our model discussed in ref. 11, and performing the integrals that remain over the fields 
<fib, <fif, an d A^, we ultimately arrive at the following "roughening model" representation 
for this partition function: 



■DC 




(12) 



n= — oo 



z= ex p 



{n b (r)} {n f (r)} 
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2 J>,(r> ±n f (r)} 2 \, (13) 



2 

r 



where n^(r) and n/(r) are integer fields that lie on the dual lattice r. 19 Also, after ap- 
plication of the Poisson summation formula, this expression may be transformed further 
into 



/oo poo ( i 

U r de b (r) U r de f (r) £ £ exp -— ^(A^ 6 ) 2 

{<fc(r)} {*/(»•)} I ^ r, M 

" 2F E( A ^/) 2 - ^ {0b + 0ff + 2nt ^ {qb9b + qfdf) }' (14) 



where <z&(r) and <?/(r) are the respective dual lattice vortex "charges" of each specie that 
range over the integers. After integrating out the fields, 6>&(r) and 6>/(r), that result from 
the Poisson summation formula, we find that the present partition function is equal to 
the product of a trivial gaussian factor with the following generalized 2D Coulomb gas 
ensemble average: 



^Coulomb = ^2 eXP 

{««(»■)} 



-(27r)^i^ ?Q (r)G^(r,r')^(0 



2 

r,r 



(15) 



where 



G aP (r,r') = (2n)- 2 [ d 2 k ^^G a p{%) (16) 

</BZ 

is the Greens function, such that the matrix inverse of its Fourier transform, G(k), is 

G ^(k) = ^-~k 2 8 a p + pq 2 g 2 u aUl3 . (17) 

The above indices a and (3 represent the internal spin/charge label, b or /, while k 2 = 
4 — 2cos k x — 2cos k y are the eigenvalues of the lattice Laplacian operator. In addition, 
u = (1, 1) represents the spin/charge component of the statistical degree of freedom [see 
Eq. (25)]. After inverting the righthand side of Eq. (17), we find that the matrix Greens 
function is decomposable into long-range and short-range components 

G aP (k) = G^(k)+G%\k) : (18) 
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such that 

where v = (1, —1) represents the spin/charge components of the physical electronic degree 
of freedom (8), and where 

\; t 2 = (qg) 2 (P b + Pf) (20) 

yields the characteristic length scale, A st , for fluctuations of the statistical gauge-field. 
The conjunction of these last few expressions (Eqs. 15-20) are, in effect, the Coulomb gas 
representation of the AH 2 model. 

The above formulae may be reduced further, however. After substituting in the long- 
range and short range components of the matrix Greens function discussed above into (15), 
algebraic manipulation then yields that this generalized Coulomb gas partition function 
may be expressed as 

^Coulomb = Yl E ex Pi 2yr ^ E - ')?ei(r)<Zei(0 + ^sr{f - f ')q' st (r)q' st (/)] 1, 

{Qb(r)}{q f (r)} I (r,r') J 

(21) 

where the physical electronic (el) flux-charge and the modified statistical (st) flux-charge 
are given respectively by 

Qei(r)=qf(r)-q b (r), (22) 
£M=(|) »M, (23) 

and where (r, r') denote combinations of points covering the dual lattice. 25 Corresponding 
to these charges are long-range (lr) and short-range (sr) potentials given respectively by 

rir(r)=/ ^(l-e^)l (24a) 

r sr (r1 = / ^-^I^tW ^ 

Jbz 2n k 2 + A s t 

Note that for the sake of technical conveniance, we have added an overall constant to the 
matrix Greens function represented by the first term in the above two equations. Due 
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to the existence of a long range force, charge nuetrality is enforced for each species and 
this constant term therefore sums to zero for all of the relevant charge configurations in 
the energy functional of the enseble average [see Eq. (15)]. Since the effective Coulomb 
gas inverse temperature scale is given by j3, we expect there to exist a binding-unbinding 
transition of the BKT type modified by short-range interactions for the physical electronic 
flux-charges, q e i(f), at a corresponding critical temperature. On the other hand, for the 
case of purely statistical flux-charge configurations, q e \(r) = and q s t(r) = <?/(r) + Q&(r) 7^ 
0, Eq. (21) indicates that there exists only short-range forces, where A st determines the 
length scale for such forces in units of the lattice spacing, a. Thus, we expect no BKT 
transition in such a case. Below, we elaborate on the consequences of the latter statements 
to show that (i) a 2D superconducting transition of the BKT type exists only for the 
physical electronic phase (8), and that (ii) all other phase auto-correlations, including 
those among the statistical phase, 



= 0/ + 06, (25) 



are short-range at all temperatures. Note that exactly the opposite behavior is exhibited 
by two square-lattice XY models stacked up on top of each other with nearest-neighbor 
coupling, where the phase difference between layers has short-range order at zero temper- 
ature, while the sum of the phases between layers shows quasi long-range order at low 
temperature. 26 

B. Auto- correlation Functions. We now wish to probe the phase correlations of the 
presently considered AH 2 model for spin-charge separated superconductitvity in two di- 
mensions. The gauge-invariant phase correlation function may be expressed as 



C12 = ^exp|z^p Q </> Q (fi) - Q (r 2 ) + q A^dr^ 



f, (26) 



(a = f, b) where the partition function Z' differs from Z only by the addition of the 
exponent — i J2 a Pa[4>a{ri) ~~ 4>a if 2) + q A^dr^] to the energy functional (4). Again, 
following refs. 18 and 19, Z' may be computed via the Villain duality transformation. 
In particular, substitution of the mathematical identity (12) that is valid in the low- 
temprature limit ultimately yields the following "roughening model" representation for 
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Z '= zZ E ex Pi -^-I][ A M^(0+P6^(0] 2 -^-I][A M n / (r)+p / r 7M (r)] i 



" J>,fr) • ///(r);n. (27) 



2 



r 



where the dual lattice "dipole" vector, ^(r), has a value of ±1 if it intersects a fixed path 
on the original lattice connecting points r\ and r*2, and it vanishes otherwize. As before, 
after the application of the Poisson summation formula, and subsequently integrating over 
the fields that this formula generates, we find that 



Z' =exp 



- 1 - p l + p -L)Ynl(r) 



r jf j, 



X 



X 



exp(-(27r)^^[ ?a (r) + ^(r)]G a Mr,r')[ ?/3 (r , ) + ^(^)]j, (28) 

{q a (r» I r,r' ) 

where the external vortex "charges" corresponding to each species are given by 

with 7](r) = A At ?7 M (r). 18 ' 19 This expression constitutes the Coulomb gas representation of 
the modified partition function Z' . 

The above expression for the modified partition function can be reduced further by 
substituting in the explicit form of the matrix Greens function G [see Eqs. (16-19)]. We 
then find that the preceeding correlation function, Z'/Z, may be factorized into 

Cl2 = Cwave(ri - ^Gcoulomb^l - f 2 ), (30) 

where 



- f 2 ) = expi - U P l + P l ) 



^v(r)[p'i^\r(r-f')+p 2 st T ST (f-f')]r](r')\ (31) 



^ PbPf fp 
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and 



Gcouiomb^i - r 2 ) = ( exp<| z ^J^ 1/2 Y^\Pe\Qei(r)T lr (f - f')+ 

+ Pst?st( r ) F sr(^- f')]l](r') 



(32) 



Coulomb 



with = ^[(/3b//3/) 1//2 p/ — (/3///9&) 1 ^ 2 Pfo] and p st = \{pf + Pb)- After some manipulation 
(see Appendix A), it can be shown that in the limit \r\ ^> A st , 



G , wave (f) = exp 



(r + 7T ^nAgt) 



exp 



v'Lt rir(r) 

2/3 7T 



(33) 



with 



Pel 



Ps't 



2/5 / 5 / 1 : 



2/3 



Pst, 



(34a) 
(346) 



and where p e i = ^(P/ — P&) an d lnA st = (2n)~ 1 J BZ d 2 k(k 2 + A^ 2 ) -1 . Hence, all phase 
correlations are short-range at all temperatures unless p s t = 0. In particular, purely 
statistical phase correlations, with p e \ = and p s t 7^ 0, are short-range. This implies that, 
within the present model, transverse statistical gauge-field fluctuations do not acquire a 
gap via a Higgs mechanism. 

A short-range order to algebraic long-range order transition of the BKT type does, 
however, exist for purely electronic phase correlations, p st = and p e \ 7^ 0, as we demon- 
strate below. But before considering this case in particular, let us first compute Gcouiomb 
for the presently considered AH 2 model on the square-lattice in general. This function 
may be calculated by following well-known methods that have been successfully employed 
in the corresponding computation for the case of the XY model. 18 ' 19 In particular, given 
the definitions 



r' 

a sr(r) = ^r sr (r - f")?7(r / ) ! 



(35a) 
(356) 
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we can re-express the latter correlation function (32) as 

Wf) 1/2 



C C ouiomb(ri -f 2 ) = ( exp<| i 2 /2 ^{\p' el q el (r)a h (r) + p st q' st (r)a sr (r)}) 



Coulomb 

(36) 



Since charge correlations in the present model are strong as a result of the existence long- 
range forces, we can make the following cummulant expansion: 

G C ouiomb(ri - r 2 ) = expi -\-^q- ([Pei9ei(r' )a h (r') + Pstq' st ( r ')vsA r ')] 

(. * r',r" 

x ^eiMairM+^M^CO])!. (37) 

But since only the lowest energy vortex-antivortex excitations need be accounted for at 
the low temperatures presently considered, the above expression then reduces to 



(38) 



Ccouiomb(^i - r 2 ) = exp 
where 

a(r') = sgn (ft - p f ) p> r (r') + min [(P b /P f ) 1/2 , (Pf/P b ) 1/2 } I&*«(r% (39) 

and where the Coulomb charge correlation function for the species labeled by i, such that 
Pi = min (p b , p f ), is given by 



<?i(0)ft(r)) = -exp<^ -2nP 



r lr (r)+minl |> & |r„(r) 



(40) 



Note, therefore, that vortex-antivortex excitations are the least energetically costly for the 
latter species. 

In the particular case of the physical electronic phase correlator, for which p st = 0, 
the above expression (38) has precisely the same form as that corresponding to the XY 
model; 18 ' 19 i.e., cr(r') = ±p e \air(r') . Due to the fact that vortices are bound to anti-vortices 
in the present low-temperature limit, we may expand cr\ r (r') and u\ r {r") with reference to 
the center of mass coordinates; i.e., 

<Tir(r') =a h (R) + ^f-Va lr (R), 
<rir(r") =a lr (R)-^f-Va lr (R), 
16 



where R is the center of mass coordinate and r is the relative coordinate between points 
r' and r". Due to the charge neutrality condition effectively imposed by the presence of 
long-range interactions, only the last terms above remain in sum over coordinates found 
in the exponent of expression (38). This ultimately results in 



Gcouiomb(ri - r 2 ) = exp 



R 



(41) 



where the factor of 1/2 found in the sum over the relative coordinate r is due to an angle 
average. However, it can be shown in a straight forward manner that 18 ' 19 



1 - 



R 



Inserting this identity into (41), we arrive at the final closed form expression 



Gcouiomb(n - r 2 ) = exp 



(42) 



Last, factoring in above the previous result (33) for the "spin-wave" contribution, we find 
that in the limit \ f\— r^l ^> A s t, the auto-correlation function (26) for the physical electronic 
phase is given by 



e iPei[<Mi)-<M2)] } ^ exp 



27T/3eff 



ri r (ri - r 2 ) 



(43) 



in the low-temperature limit, with an effective temperature scale 

^*0>)«(r)>. 



(44) 



\r\>a 



We now recall that T\ T (r) = ln(2 3 / 2 e 7 r), where 7 is Euler's constant. 19 Hence, in the limit 
A s t — >• 00, where statistical gauge-field fluctuations in the energy functional (4) are "frozen 
out" , inspection of Eq. (40) reveals that j3~^ diverges at 2n min(/?&, (3f) = 4. This implies a 
BKT transition temperature of T c = ^ min( J5, Jf) below which algebraic long-range order 
sets in (see dashed lines, Fig. 1). The latter result is precisely what one expects from 
employing the Ioffe-Larkin formula for the London penetration length in a spin-charge 
separated superconductor; 3 ' 8 i.e., the mean-field approximation result (2). On the other 
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hand, in the general case of finite A st , similar considerations reveal that (3~ s diverges at 
2tt/3 = 4, which implies a BKT transition temperature given by Eq. (9). Hence, both the 
previous continuum limit analysis and the present low-temperature analysis arrive at the 
same answer concerning the nature of the BKT phase-transition in the AH 2 model. 

Last, we note that the previous analysis for the Coulomb gas factor, Gcouiomb, in 
the particular case of the physical electronic phase correlator (p st = 0) is also valid for 
the general case in the limit of maximum statistical gauge-field fluctuations, A st — ■> [see 
Eq. (32)], 18 ' 19 with the expception that we must replace p e \ in Eq. (42) with the orginal 
expression (34a) for p" v Factoring in the latter result with the previous result (33) for the 
"spin-wave" contribution, one finds that the correlation function (26) is given by 



Hence, we see that in the limit of maximum fluctuations in the statistical gauge-field 
considered here, while the phase auto-correlations are generally short-range , a remnant of 
the long-range behavoir that exists in the special case of the electronic auto-correlations 
function (43) persists. 



We have demonstrated above that the statistical gauge-field does not acquire a gap via 
the Higgs mechanism, since the corresponding auto-correlations for the statistical phase 
(25) are short-range. However, a gap in a U (1) gauge-field can also result from confinement 
effects in 2+1 dimension, 27 for example. Such effects are conveniantly probed by the Wilson 



a perimeter law for this quantity at all temperatures below the BKT-transition discussed 
above, which simply reflects the existence of bound vortices, and which corresponds to a 
"deconfined" phase. Also, although we do find that the Wilson loop exhibits a "confining" 
area law at temperatures above the BKT phase-transition, we argue that this behavior is a 
general property of pure Abelian gauge theories in two euclidean dimensions. We thereby 
conclude that the statistical gauge-field excitations remain gapless at all temperatures in 
the present AH 2 model. 




(45) 



V. Statistical Gauge-field Excitations and the Wilson Loop 



loop, (exp(ze' § A^dx^)). Below we show that the presently considered AH 2 model displays 
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A. Low- Temperature Phase. Let us first consider the Wilson loop for a large contour C 
in low temperature limit. As in the previous case of the phase auto-correlation functions, 
it can be expressed as 

W(C) = ^exp ^e' £ A^dr^i ^ = J' (46) 

where the partition function Z' differs from Z only by the addition of the exponent 
—ie'§ c A^dr^ to the energy functional (4). Again, following refs. 18 and 19, Z' may 
be computed via the Villain duality transformation by substitution of the mathematical 
identity (12), which is valid in the low-temperature limit. Generalizing the treatment of 
Jones et al., 11 we find ultimately that Z' has the following "roughening model" represen- 
tation: 

E E exp{-^-E[ A M^(r)] 2 -^-Et A ^/W] 2 

{n b (r)} {n f (r)} I Ph r,n W r,n 

- E h (r) + n/(r) + Vq J(r) ] 2 r (47) 

r yy J 

where J(r) has a value of 1 if the point r is within the contour C, and it vanishes otherwize. 
Again, after the application of the Poisson summation formula, and subsequently integrat- 
ing over the fields that this formula generates, we find that Z' is given by the product a 
trivial "spin-wave" factor with the following Coulomb gas enesemble average: 



^Coulomb = E eXP ] -(^^E 

{<la(r)} I r,r' 



q a (r)G al3 (r,r')qf3{r') + in 1 e'qgq ce (r)G i ^J (r,r')upJ(r') 



(2n)- 2 e' 2 (qg) 2 J(r)u a G% ) (r,r')u p J(r') + (2n)~ V 2 ^" 1 J(r)5 r y 



(48) 



again where u = (1, 1) represents the spin/charge components of the statistical degree of 
freedom (25), and where G^(r, r') represents the short-range component of the matrix 
Greens function [see Eqs. (16-19)]. The above expression constitutes the Coulomb gas 
representation of the modified partition function Z' . 

To further reduce expression (48) for ^coulomb' we recan that 
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Hence, in the limit of large contours, C, the sum over the relative coordinates in the third 
term of this expression can be extended to the entire lattice. In this case we find that the 
last two terms in (48) sum to zero; i.e., 



2„/2 



(27r)- 2 e 



_ E |g (A+ ^ J(r)+E a + _L )j(r) 



by (20). Note that the first term above corresponds to the third term in (48), where the 
sum over the relative coordinates has been already performed. Similarly, since (r, r') 
is short-ranged, we have that in the limit of large contours, 

J2q«(r)G ( :$(r,r')upJ(r') = ^f 1 ^ £[&<&(r) + /3 f q f (r)] J(r). 

r,r' PbPf r 

The lefthand side of this expression is proportional to the second term in (48), which 
implies that the Wilson loop, Z'/Z, is given by 

W(C) = /exp -2Kie'qg\lY,{(3bq h {r)+(3fqf{r)]J{r) \ . (49) 

\ L r \ I Coulomb 

Since vortices are bound to anti- vortices at low temperature due to the existence of long- 
range interactions, we may employ a cummulant expansion to compute this average, as in 
the previous calculation of phase correlators [see Eqs. (36) and (37)]. Such an approxima- 
tion yields 

W(C) = exp l-(2Tve'qg) 2 X 4 st \ ^ ( [(3 b q b (r') + (3 f q f (/)] [(3 b q b (r") + (3 f q f (r")]> J (/) J (r") 



= exp< 



2ne' 

qg 



1 + max| ^'^ 



-]>>(rW0>J(rV(O , (50) 



again where the species label i is such that Pi = min(,9 fc , f3f). Since the charge correlations 
decrease rapidly, we may expand J{r') and J(r") with reference to the center of mass 
coordinates; i.e., 

J(r') =J(R) + l?.VJ(R), 
J(r") =J(R)-±f-VJ(R), 
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where R is the center of mass coordinate and r is the relative coordinate between points 
r' and r". Due to the charge neutrality condition effectively imposed by the presence of 
long-range interactions, only the last terms above remain in sum over coordinates found 
in the exponent of expression (50). This ultimately results in 



inn ex,^'' 



exp< 




, ft Pf 

1 + max| /Vft 



1 + max I 



A, §i 
P f ' A, 



E^(«(°)«( r ))Ez^ J ( fl )] : 

r 

5>%((%(r)> 



P 



(51) 



where that prefactor of \ in the sum over the relative coordinates arises from an angular 
average. Above, P denotes the length of the contour C. The latter expression represents 
the final result of our manipulations. 

We observe, therefore, that at the low temperatures presently considered, the Wilson 
loop exhibits a perimeter law. And since the vortex charge correlator, (qi(0)qi(r)) , corre- 
sponding to the species with the least costly vortex excitations is given by Eq. (40), we 
also observe that the coefficient to the above perimeter-law diverges precisely at the BKT 
phase transition temperature (9). This suggests that a "confining" area law behavior exists 
above T c , which we in fact argue for below. Last, it is also of interest to remark that this 
coefficient vanishes in the limit Jb 3> J/ or J/ 3> Jb- 

B. High- Temperature Phase. Clearly, in the high-temperature limit, /?& = and 
Pf = 0, the present AH 2 reduces to a pure Abelian gauge theory, which is trivially "con- 
fining" in 1+1 dimensions, and which therefore shows an area law for the Wilson loop. 
Also, in either the limits /3b = or (3f = that correspond to the one-component model, 
it is well know that the Wilson loop follows an area law at all temperatures in the case 
of fractional probe charges e'. 11 It has recently been argued, however, that such "con- 
finement" is trivially due to pure 1+1 dimensional statistical "electro-magnetism" that 
is renormalized "dielectrically" . 21 In this case, one finds that gauge-field fluctuations are 
characterized by a free U(l) action [see the last term in Eq. (4)], with a renormalized 

1/2 

"charge" , g/e s l , that is on the order of g for temperatures above a cross-over temperature 
scale, 8 T„ 



l c/o ~ Js (s = 6 or/), and that becomes exponentially small below this cross-over 
temperature scale. 21 This ultimately leads to a cross-over phenomenon between a "strange" 
metal phase at high temperatures and its absence at low temperature. 8 Given that the only 
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true phase transition in the present AH 2 model is the BKT transition itself, by continuity 
in the phase diagram (see Fig. 1), (i) the Wilson loop should follow a "confining" area 
law at all temperatures above T c for fractional probe charges e' and (ii) we expect that a 
similar crossover phenomon occurs in the present AH 2 model in this temperature regime, 
with the exception that the crossover temperature scale is given by T c / Q ~ max( J/) in 
this case. 

In summary, the Wilson loop (46) of the statistical gauge-field in the presently con- 
sidered AH 2 model shows a perimeter law at low temperatures T < T c , such that the 
corresponding coefficient diverges at T c . Also, in the general case of fractional probe 
charges, e' ^ nq (n = 0, ±1, ±2, ...), the Wilson loop shows an area law for tempera- 
tures above the BKT transition, which we argue implies a gapless spectrum for statistical 
gauge-field excitations. Since we showed in the previous section that there exists no Higgs 
mechanism for the statistical gauge-field to acquire a gap due to the complete absence of 
long-range order in the statistical phase (25), and since the perimeter law exhibited by the 
Wilson loop in the low temperature phase does not indicate the existence of short-range 
"electromagnetic" interactions, we argue that the statistical gauge-field remains gapless as 
well in this phase. The occurence of a perimeter law, therefore, is simply a signal of the 
vortex binding-unbinding transition present in the system. 



VI. Universality Class 



By consideration of both the continuum limit (section III) and the low-temperature 
limit (section IV), we have shown above that a short-range order to algebraic long-range 
order phase transition occurs at a transition temperature given by Eq. (9) for the physical 
electronic phase (8) in the present AH 2 model for spin-charge separated superconductivity. 
But what is the nature of this transition in physical terms? Below, we argue that the fore- 
mentioned transition is in the same universality class as that of the XY model. Hence, we 
expect a corresponding universal jump in the superfluid density. 16 

Let us now derive the renormalization group equations for the BKT transition in the 
physical electronic phase of the present model by closely following the treatment of the 
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corresponding problem in the case of the XY model. 19 From eq. (44), we have that 

/ \ 2n(3(r)-3 

1 1 o 2 f°° dr ( a 



&ff (3 



a \ r 



(52) 



where by Eq. (40) for the charge correlations, (3{r) is a smooth function satisfying 



A; r < A 



St 



/3(r) 



/?; r > A 



St 



(53a) 
(536) 



and where the activity at r = r is given by 



y = exp 



a 



-7r/?(r ) In— 



(54) 



Let x = r/a, and suppose that we make the following approximation to the integral 
appearing in the rigthhand side of Eq. (52): 

"A s t 



i; 



dxx 



3-27T/3; 



+ / dxx 



3-27T (3 



A 



4-27r/3i 
st 



- 1 



POO 

In A st + [1 + (4 - 2tt/?)1ii A st ] / dxx 3 ' 2 ^ 



in the limit A s t — > 1. Hence, in this limit, we have by (52) that 
1 1 



Pefi 



= + 7r3 y 2 ln Agt + ^3^2 [i + (4 _ 27r/ g)ln Ast] / dxa; 



3-2tt/3 



This formula may be rexpressed as 

1 1 



where 



; + ^y\ / ete 3 " 2 ^, 

Peff P\ Jl 



1 1 3 2, x 

Px P 

y x = y+(2- 7r/3A)ylnA st . 
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(55) 



(56) 



(57) 



(58) 
(59) 



The differential form of the latter two relations are simply the Kosterlitz renormalization 
group equations 16 ' 19 

KtJ-K^^vl (60) 



dX 



St 



\ s t^yx = (2-n[3 x )y x , (61) 

with the length renormalization scale determined by the model parameter A s t- 

In summary, given the validity of approximation (55), we find the BKT transition 
obtained previously in the present AH 2 model for the physical electronic phase is in the 
same universality class as that of the XY model. We therefore expect a universal jump in 
the superfluid density, following that predicted for the XY model. 16 



VII. Summary and Discussion 

In this paper, we have introduced a square lattice gauge theory model for spin-charge 
separated superconductivity in the presence of statistical gauge-field excitations. Such 
excitations appear naturally in spin-charge separated descriptions of strongly interacting 
electrons systems, and they represent so-called chiral spin fluctuations. 10 After analysing 
this model in various limits (continuum and low-temperature), we found that it explicitely 
shows a superconducting transition of the BKT type for the physical electronic phase (8) 
at a critical temperature given by (9). As expected, we see that statistical gauge field 
excitations suppress the superfluid transition temperature with respect to mean-field; i.e., 
T C <T C (0) . It was also argued in the preceeding section that this transition shares the same 
universality class as that of the XY model, and hence that we expect a corresponding uni- 
versal jump in the superfluid density. All other phases, however, were found to show only 
short-range auto-correlations at low-temperature. In particular, the fact that statistical 
phase (25) correlations are short-range at all temperatures implies that the corresponding 
transverse statistical gauge-field fluctuations do not acquire a gap via the Higgs mecha- 
nism. Furthermore, we have argued that because of dimensionality, neither "confinement" 
effects produce a gap in these excitations. Hence, we claim that the statistical gauge-field 
excitations remain gapless at all temperature in the present model. It is important to note, 
however, that gapless statistical gauge-field excitations also have been shown to exist at 
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the meanfield level in the commensurate flux phase of the t — J model, which is thought 
to be a spin-charge separated anyonic superconductor, but they were to found acquire a 
gap once dynamical effects were included. 22 It is therefore not certain that such gapless 
excitations will persist in the present model once such effects are accounted for. We hope 
to address this issue in a future publication. 

The superfluid transition found in the AH 2 model for spin-charge separated supercon- 
ductivity can be understood a posteriori as follows. Let ^/ = and ^i, = be 
the respective order parameters of each specie. Now suppose that the superconductivity 
in this spin-charge separated system results from conventional Cooper pairing of electrons, 
with an order parameter given by = (cjf&]cj' ib\,) . Then the mean-field approximation 
yields that the true superconducting order parameter satisfies 

It is thus not surprising to find that the physical electronic phase (8) develops quasi-long- 
range order in the present model. What is indeed important to note, however, is that 
fluctuations in the statistical gauge field do not suppress the transition altogether. 

Yet can the present model explain, or at least describe, the phenomenon of high- 
temperature superconductivity? Mean-field RVB treatments of the t — J model yield 
a superfluid phase stiffness for spinon pairs of J/ ~ J(l — x/xq), where x denotes the 
hole concentration and xo denotes the critical concentration beyond which there is no 
pairing instability. 9 Also, since we have assumed throughout that the holons form pairs, 
then their corresponding superfluid phase stiffness is given by Jb ~ t'x, where t' denotes 
the hopping matrix element for such pairs. Now suppose that the latter hopping matrix 
element satisfies t' ~ 1000 K, which is conceivable since the corresponding matrix elements 
for a single electron is on the order of t ~ 5000 K. Then since Eq. (9) implies that the 
superconducting temperature satisfies T c ~ mm(t'xo, J), and since J ~ 1000 K while the 
maximum hole concentrations for oxide superconductivity are typically near xq = 0.2, we 
have that T c ~ t'xo ~ 200 K, in agreement with experiment. 1 More importantly, however, 
Eq. (9) for the critical temperature coupled with the above functional dependences for 
the superfluid stiffnesses of each species with hole concentration imply a phase-diagram for 
the superconducting versus normal phase that is shown in Fig. 1. The shape of the phase 
boundary qualitatively resembles that of the high-temperature superconductors, 20 but it 
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differs substantially in shape from that predicted by Ginzburg-Landau treatments of spin- 
charge separated superconductivity. 8 Finally, since g 2 = T/xd in this case, where Xd = 
Xf + Xb ~ J(l — x) +tTb/T is the sum of the diamagnetic susceptibilities of each species, 5 ' 6 
with Tb ~ tx, we then have that A s t = q~ x [Xd/iJb + Jf)] 1 ^ 2 ~ 1- Hence, the renormalization 
group analysis discussed in the previous section is valid, and we therefore expect a universal 
jump in the superfluid density at the critical temperature in the absence of interplanar 
coupling. Note that the present description of oxide superconductivity differs substantially 
from that proposed by Anderson and coworkers, 28 where it is assumed that interlayer 
interactions drive the phenomenon. Although it is proposed here, on the contrary, that 
oxide superconductivity is primarily a one-layer effect, we nevertheless hope to study the 
problem of a few coupled layers of square lattice AH 2 models in a future publication. It 
is quite possible that the nature of the presently discovered superfluid transition changes 
dramatically once the model is extended into the third dimension, which could be relevant 
to the layered high-temperature superconductors. 8 

Finally, what is the nature of the normal state above T c ? It was previously remarked 
in section II that separate "strange" metal, Fermi-liquid, superconducting, and spin-gap 
phases exist in certain spin-charge separated treatments of the t — J model in two dimen- 
sions at the mean-field level, 8 as shown by the dashed lines in Fig. 1. We remind the reader 
that the "strange" metal phase corresponds to the absence of superfluidity in both species, 
the Fermi-liquid phase to the appearence of superfluidity in the holon species alone, the 
superconducting phase to the appearance of superfluidity in both species, and the spin- 
gap phase to the existence of superfuidity in the spinon species alone. In this paper, we 
have demonstrated explicitely that the only true phase transition that remains beyond the 
meanfield approximation is the superconduting one depicted by the solid line in Fig. 1. In 
particular, fluctuations in the statistical gauge field representing chiral spin-fluctuations 10 
destroy the meanfield transitions between the "strange" metal phase and the Fermi-liquid 
phase, and between the "strange" metal phase and the spin-gap phase. 8 This effect can be 
studied more easily by considering the case when only one of the species has a superfluid in- 
stability; i.e., the one-component Abelian-Higgs model corresponding to /?& = or f3f = 0. 
It has been explicitely shown in this case, by using techniques similar to those employed in 
this paper, that the BKT superfluid transition present in the absence of statistical gauge- 
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field fluctuations becomes a cross-over phenomenon. 21 The cross-over temperature scale 
satisfies T c / Q ~ T s , where T s denotes the meanfield transition temperature of the species 
in question, and it vanishes in the limit of maximum gauge-field fluctuations, A st — > 0. By 
continuity with the present AH 2 model, we expect that this crossover phenomenon persists 
in the normal phase above T c ; i.e., T c / Q ~ max(Tj, Tf). Furthermore, in the limit A s t — > 0, 
the normal phase should be entirely "strange" , with a characteristic T-linear resistivity. As 
discussed above, however, reasonable estimates based on the t — J model yield a statistical 
length scale satisfying A st ~ 1. Hence, we expect that the cross-over phenomenon from 
the "strange" metal regime into the Fermi-liquid and spin-gap regimes should persist in 
strongly interacting 2D electron systems. Such a cross-over scale (see the upper lines in 
Fig. 1) is consistent with the narrow doping window for linear-in-T resistivity recently 
observed in the superconducting oxides. 29 
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Appendix A: "Spin- Wave" Component of Phase Correlator 



To calculate the "spin-wave" component of the phase auto-correlation function (31) 
in the AH 2 model at low-temperature, given by 

G wave (ri2) = expX wave , (A.l) 



with 



A wa ve 

1 P 



Y,v(r){p?ir h (f-f')+pir si (f-f')] V (r'), (A.2) 



we will follow the standard treatment of the corresponding problem for the XY model. 18 ' 19 
Let us suppose that r*i 2 = r 1 — r 2 is directed along the x-axis. Then clearly, 

---Jtt E Pd[2ri r (x-x',o)-ri r (x-x , ,i)-r lr (x-x , ,-i)] 
1 3 ri2 ~ 1 

--1TT E Pl[^sr(x-x\0)-T ST (x-x\l)-T sr (x-x\-l)]. (A3) 

But since T\ r (r) and r sr (r) — r sr (oo) are proportional to the respective long-range and 
short-range lattice Greens functions, then by definition [see Eqs. (16), (18) and (19)] they 
satisfy 

4T lr (x, y) - [T h (x + l,y)+ T h (x - 1, y) 

+T h (x, y + l) + r b (x, y-l)} = -2TrS Xt0 6y,o, (A A) 



4r sr (x, y) - [T SI (x + l,y) + T SI (x - 1, y) 

+r sr (x, y+l) + T sr (x, y - 1)] + A" t 2 [r sr (x, y) - r sr (oo)] = -2^,0^,0- (A.5) 

Substituting in the latter identities into (A. 3) for the Greens function terms at points off 
of the y = axis, and then using the fact that the resulting sum over x and x' is a partially 
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telescoping series, 18 ' 19 we find that 
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2(3 /2 r lr (r 12 ) 

"Pel" 
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ri2 — 1 



2/3 2 I r sr (n 2 ) , x -2 1 r-p / n -p / /\] 
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But in the limit r\ 2 ^> A st , we have the identity 

ri2 — 1 -i r 12 — l oo 



A s"t 2 E ^[r sr (oo) - r sr (o; - x')} \~ t 2 E E^[ r -(°°)- r -( a; - a;/ )]=ri2 



x,x'=0 



x=0 x'=0 



Substituting this relationship into the last term of (A. 6), we find by (A.f) that the "spin- 
wave" component of the phase correlator in the AH 2 model on the square-lattice is given 
by 



G W ave(^12) = exp 



2(3 



PbPf 



Pst( r 12 + 7T /nA st ) 



exp 



2(5 /2 r lr (r 12 ) 



(A.7) 



in the limit |f| ^> A s t, where 



/nA st = r sr (oo) = (27T)" 1 / d 2 k(k 2 + A" 2 )" 1 . 

Jbz 



(A.8) 



Note, of course, that the function defined above coincides with the Napierian logarithm in 
the limit A st ^> 1. 

Appendix B: Phase Correlations in One-component Model 

Consider the gauge-invariant phase correlation function 

C 12 = (exp{ip 6 [0 6 (n) - 4> b {r 2 ) + q [ 2 A^dr^]}) = Z'/Z (B.l) 

Jri 

in the case of the one-component Abelian-Higgs model (3f = 0, where the partition function 
Z' differs from Z only by the addition of the exponent —ipb[4>b{ r i) ~ 4>b{ r 2) + Q A^dr^ 

29 



to the energy functional (4). Then by (27), we have that the Villain duality transformation 
of this partion function is given by 

z ' = E ex p{-^-Et A M^w+^w] 2 -^E^( r )}- ( s - 2 ) 

{n b (r)} b r,n r 

As discussed in section IV, n&(r) ranges over the integers and r covers the dual lattice, 
while the dual lattice "dipole" vector, ^(r), has a value of ±1 if it intersects a fixed path on 
the original lattice connecting points r\ and r*2, and it vanishes otherwize. 18 ' 19 Therefore, 
in the limit A st <C 1, where statistical gauge- field fluctuations are maximized, the trivial 
configuration n^r) = dominates the ensemble average above, resulting in 

Cw^exp^^lfi-fal). (fl.3) 

Notice, as expected [see Eq. (9)], that this correlation function shows no evidence for a 
phase-transition at non-zero temperature. 

It was also shown previously that the application of the Poisson summation formula 
to the "roughening model" representation expression (B.2) leads ultimately to the factor- 
ization C12 = C wave (ri — r^Gcouiombtji — t 2 ) of the gauge-invariant phase correlation 
function. In the present case of (3f = and pf = 0, Eq. (A. 7) yields that 

2 

Cwave^ -f 2 ) ^expf-^-df! -f 2 | +7r- 1 ZnA st )l (B .4) 

for r 3> A st , while Eq. (36) yields that 

Gcouiomb(ri - f 2 ) = (exp\ip b q b (r)a ST (r) }) , (B.b) 

\ >- ^— -" ' J I Coulomb 

r 

where the corresponding partition function (21) for the screened Coulomb gas ensemble is 
given by 25 

■^Coulomb = ^ ] © X P 

{ 9i) (r)} I (r,r') J 

Unlike the previous treatment in the case of the AH 2 model, however, we now take the 
original definition 



d h ^ 1 

. - 2 



'BZ ^ /C 2 + X st 
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for the interaction energy, since charge conservation is no longer guaranteed due to the 
absence of long-range forces. Following ref. 11, in the limit A st 3> 1, in which case the 
interaction between vortices is negligible at large distances, we may approximate the above 
enemble average by 



^Coulomb — II r ^2 eX P 

Qb(r)=-oo 



-2irf3 b ^(ln\ st )q 2 b (r) 



(B.8) 



Hence, if we truncate the sum over vortex charge to qb{r) = 0, ±1, then by (B.5), 

^ IT r {l + 2exp(-7r/jb/nA st )cos [p b cr S r(r)}} 
Coul ° mb " n r [l + 2exp(-7r/3 6 /nA st )] 

^ exp(2e-^ bZnAst ^{cos [p b a sr (r)} - 1}) 

r 

exp(2\-^ b ]T{cos [p b a sr (r)] ~ 1}) 

r 

-> 1 

in the limit A st — > oo and (3 b — > oo. Thus, in these limits, the gauge-invariant phase 
correlation function, C12, is approximately given by the "spin- wave" contribution (B.4) 
for r > A st . 

In the low-temperature limit of the one-component Abelian-Higgs model, therefore, 
we generally find that the coherence length for the holon phase correlations is given by 
£,b ~ Pb = Jb/T in units of the lattice constant, a, contrary to previous claims in the 
literature of an exponentially divergent length. 8 
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Figure Caption 



1. Shown is the phase boundary between the superconducting and the normal phase as a 
function of hole carrier concentration, x. The superconducting transition temperature 
is determined by Eq. (9), with f J b = (200 K) • x/x and § J f = (800 K) • (1 - 
x/xq), where xq ~ 1 denotes the critical concentration above which the spinon pairing 
instability is absent. The upper set of dashed lines represent cross-over regions below 
which "strange" metallicity begins to dissappear. 
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